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Abstract 



A relativistic mean-field model of nuclear matter with arbitrary proton frac- 
tion is studied at finite temperature. An analysis is performed of the liquid- 
gas phase transition in a system with two conserved charges (baryon number 
^ and isospin) using the stability conditions on the free energy, the conservation 

laws, and Gibbs' criteria for phase equilibrium. For a binary system with two 
^-H phases, the coexistence surface (binodal) is two-dimensional. The Maxwell 

construction through the phase-separation region is discussed, and it is shown 
that the stable configuration can be determined uniquely at every density. 
Moreover, because of the greater dimensionality of the binodal surface, the 
liquid-gas phase transition is continuous (second order by Ehrenfest's defini- 
tion), rather than discontinuous (first order), as in familiar one-component 
systems. Using a mean-field equation of state calibrated to the properties of 
nuclear matter and finite nuclei, various phase-separation scenarios are con- 
"'T' sidered. The model is then applied to the liquid-gas phase transition that 

O may occur in the warm, dilute matter produced in energetic heavy-ion col- 

^ lisions. In asymmetric matter, instabilities that produce a liquid-gas phase 



separation arise from fiuctuations in the proton concentration (chemical in- 
stability), rather than from fiuctuations in the baryon density (mechanical 
instability). 



1 



I. INTRODUCTION 



The determination of the properties of nuclear matter as functions of density, tempera- 
ture, and the ratio of protons to neutrons is a fundamental problem in nuclear physics. To 
achieve this goal, one must study not only the ground and excited states of normal nuclei, 
but also highly excited nuclei created in nucleus-nucleus collisions and nuclei far from sta- 
bility, which may be created in radioactive beams. In this work we consider the properties 
of equilibrium nuclear matter at finite temperature and arbitrary proton fraction, and we 
describe some new qualitative features that may be relevant for the energetic collisions of 
heavy ions. 

There have been numerous theoretical studies of the dynamics of medium-energy heavy- 
ion collisions [1-6] . Some of these are based on equilibrium thermodynamics and focus on the 
nuclear matter phase diagram [7-11]. The basic feature is the liquid-gas phase transition at 
low densities and moderate temperatures, and how the nuclear system evolves through vari- 
ous phase-separation boundaries (binodals) and instability boundaries (spinodals) [7,12,13]. 
Other analyses concentrate on the nonequilibrium evolution of the system by describing the 
early stages of the collision through cascade simulations [14,15], for example, or the late 
stages of the collision using models of fragmentation or nucleation [16-18]. Although the 
equilibrium analysis oversimplifies this problem and is applicable only to certain aspects of 
the evolution, it is useful for providing an orientation to the complex dynamics by giving 
concrete descriptions of the phase-separation process, which may be difficult to character- 
ize in more microscopic formulations. Here we will follow the thermodynamic approach 
and focus on the qualitatively new features that arise when a system with two conserved 
charges (i.e., baryon number and electric charge) undergoes a liquid-gas phase transition. 
These features have been discussed earlier by Barranco and Buchler [19], who use a simple 
phenomenological equation of state with an interaction energy that is temperature inde- 
pendent, and by Glendenning [20], who focuses on zero-temperature matter and neutron 
stars. More recently [21], Pethick, Ravenhall, and Lorenz used a similar analysis to study 
the composition of neutron star crusts. 

The main ingredient in the analysis is the nuclear equation of state, and ours is based 
on a relativistic mean-field model involving the interaction of Dirac nucleons with scalar 
and vector mesons [22,23]. There are several reasons for choosing this type of model. First, 
recent developments in the application of these models to the structure of nuclei show that 
they can provide an excellent description of bulk nuclear properties throughout the periodic 
table, provided that nonlinear scalar and vector self-interactions are included [23-34]. In 
fact, relativistic mean-field models describe these properties as well as or better than [35] 
any other microscopic model presently available. Thus we have a way to calibrate our 
equation of state at zero temperature and normal nuclear densities, and then extrapolate 
into the warm, dilute regime appropriate for the phase transition. Second, the mean-field 
approximation is known to be thermodynamically consistent; that is, it satisfies the relevant 
thermodynamic identities and the virial theorem [36]. Moreover, although the mean-field 
approximation oversimplifies the nuclear dynamics, it allows for easy computations and is 
commensurate with our simplified description of the collision dynamics. Finally, the ability 
to determine the model parameters analytically from a specified set of zero-temperature 
nuclear matter properties [30,34] allows us to easily study variations of the model and to 
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assess the sensitivity of our results to the nuclear compressibility and symmetry energy, both 
of which are not particularly well known. 

The models we study involve nonlinear scalar and vector self-interactions through fourth 
order in these fields, as first proposed by Boguta and Bodmer [24], and later generalized 
by Bodmer and Price [28,30]. The desired nuclear symmetry energy is achieved using the 
simplest possible coupling of the p meson to the nucleon, and there are no p self-interactions. 
Although it is known that this way of generating the symmetry energy is a simplification 
of more complex mechanisms (for example, one-pion-exchange in a Hartree-Fock calcula- 
tion produces significant contributions [37]), and that the dependence of the energy on the 
isovector density [p^ = pp — p^) may be more complicated when one is far from symmetric 
matter, the qualitative features studied here arise because the symmetry energy is repul- 
sive for all values of /J3, and this repulsion grows monotonically as \p3\ increases. Thus the 
present simple description is adequate for our purposes, and we expect that other models 
of the symmetry energy with these features will produce similar results. Moreover, we can 
adjust the symmetry energy to any reasonable value, and it is straightforward to generalize 
our model, which we leave as a topic for future investigations. 

Similar relativistic mean-field calculations of warm, symmetric nuclear matter were per- 
formed earlier by one of us [36], but there the original Walecka model [38] was used, which 
allows only a rough calibration to observed bulk nuclear properties. Here we want to ex- 
plore arbitrary proton fractions and to present results for the best relativistic equations of 
state available. Moreover, this earlier work concentrated on the Lorentz covariance of the 
relativistic mean-field theory; the same analysis could be applied to the class of models 
studied here, but this is unnecessary, and we carry out all calculations in the rest frame of 
the warm matter. Indeed, even the use of a relativistic mean-field theory is not essential for 
our discussion; any equation of state that is accurately calibrated to nuclei should produce 
similar results. 

We emphasize that although we present quantitative results based on an accurate nuclear 
equation of state, our focus is on the new qualitative features that arise in liquid-gas phase 
transitions in systems with more than one conserved charge. We believe that our simplified 
hydrodynamic, thermodynamic, and mean-field description of nucleus-nucleus collisions is 
the most transparent way to do this, and we expect that signatures of these new features 
will survive in more sophisticated calculations. In particular, our primary result is that the 
liquid-gas phase transition in asymmetric nuclear matter is of second order, in contrast to the 
familiar first-order (van der Waals) transition that occurs in one-component systems (and 
in symmetric nuclear matter). If such a phase transition actually occurs in medium-energy 
heavy-ion collisions, then generically it should be smoother than a first-order phase transi- 
tion. To be more precise, even if the assumption of thermodynamic equilibrium is totally 
justified, and even if the finite-size effects on the phase separation are totally negligible, we 
show that the generic liquid-gas phase transition must he continuous] for example, it occurs 
over a temperature range of several MeV for matter that is roughly 40% protons. Thus 
observable signals of this transition should be more continuous than previously expected 
from conventional analyses, which find a discontinuous (first-order) phase transition in the 
thermodynamic limit [16]. 

Furthermore, we find that in asymmetric matter, the relevant spinodal that signals in- 
stability to infinitesimal fiuctuations arises from variations in the proton concentration at 
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fixed pressure. This is in contrast to tlie usual scenario (and tlie one relevant for symmetric 
matter), in which spinodal decomposition occurs due to fluctuations in the baryon density 
at points of mechanical instability. Thus microscopic models that attempt to describe this 
decomposition using nucleation or fragmentation must be generalized to allow for different 
proton concentrations in the liquid and gas phases. 

The outline of this paper is as follows: In Sec. II, we present a general discussion of 
the thermodynamics of phase transitions in multicomponent systems. We emphasize that 
the stability conditions, conservation laws, and Gibbs' criteria produce a set of equations 
that completely specifies the thermodynamic variables. (That is, there are equal numbers 
of equations and unknowns.) While some of this material has been presented elsewhere, we 
believe this is the first derivation that proceeds directly from a study of the free energy of 
the system. In Sec. Ill, we describe the relativistic mean-field model, present the relations 
that determine the equation of state, and discuss how the model parameters are specified. 
Some numerical procedures used to obtain our results are also discussed. Section IV deals 
with a study of the binodal surface; various phase diagrams and Maxwell constructions 
are illustrated, and the evidence for the second-order liquid-gas transition is presented. In 
Sec. V, we apply our simple model to heavy-ion collisions and describe the phase coexistence 
region and the various spinodals. We also consider how the continuous phase transition may 
affect the evolution of the system. Section VI contains a short summary. 

As a guide for the reader, we note that the material in Sees. II through IV is presented 
before the applications in Sec. V for both logical and pedagogical reasons. Nevertheless, 
because of the formal nature of the initial discussion, it may be useful on a first reading to 
consider Sees. II through IV briefiy and then to concentrate on the results of the formalism 
that are displayed in Sec. V. After gaining some insight into the behavior of a specific 
two-component system, one can then return to the formal material and consider it in more 
depth. 
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II. PHASE TRANSITIONS IN MULTICOMPONENT SYSTEMS 



We consider a system characterized by a hamiltonian H and a set of n mutually com- 
muting charges Qi. Here a conserved charge does not necessarily imply an independent 
particle species; it includes any conserved quantity resulting from an underlying symmetry 
of the system, for example, electric charge, baryon number, total angular momentum, etc. 
Moreover, particle numbers may change during a process (by decay, for example), and it is 
actually the conserved charges that are relevant for the thermodynamic analysis. 

The equilibrium state of the system enclosed in a volume V is completely described by 
the thermodynamic potential 



O(r,y,/i0 = -^lnTrexp(^-/3(i7-^/i40) , (1) 

where /3 is the inverse temperature. Thus the average charges and the pressure can be 
obtained from 

g,(r,y,/.o = -(|^) , (2) 

piT,V,^,) = -(^) . (3) 
\oVJt,{„,} 

In some regions of T and /ii, several choices of the fii may lead to the same p, and this 
allows for the possibility of phase transitions. To be more specific, let us first perform a 
Legendre transformation and consider the Helmholtz free energy F. We assume a system 
with specified charges Qi in a large volume so that surface effects can be neglected. We 
then have 

n{T,V,iJ^i) = -Vp{T,iJ^i) , (4) 
Qz = Vp^{T,fi,) , (5) 

and correspondingly, 

F{T, V, go = vj^{T, Pi) = o(r, y, pi) + y ^ p,p,{T, pi) , (6) 

i 

with 



The system will be stable against separation into two phases if the free energy of a single 
phase is lower than the free energy in all two-phase configurations. This requirement can be 
formulated as 

^(T, Pi) < (1 - A)^(r, p',) + A^(r, (/[) , (8) 

with 
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= (1 - \)p[ + \p'l , < A < 1 , 



(9) 



where the two phases are denoted by a prime and a double-prime. In formal terms, stability 
implies that the free energy density is a convex function of the densities pi [40]. Convexity 
implies that stability against separation into two phases also guarantees stability against 
separation into an arbitrary number of phases (which we denote by superscripts a, /3, . . . , 
Lo). That is, by Jensen's inequality on convex functions, Eqs. (8) and (9) are equivalent to 

HT.Pi)<Y.^»HT.p-) , (10) 

{«} 

with 



Pt = J2^"Pt = 1 ' (11) 

{«} {«} 

where the sums are over all phases {a}. The parameters A" = V" /V specify the volume 
fraction occupied by each phase. The second equation (9) or (11) ensures that the overall 
charges are conserved: 

vp. = j2yy"^ (12) 

{«} 

with V = J2V". (13) 

{«} 

Equations (8) and (10) are global criteria for the stability of the one-phase system. If 
these inequalities are satisfied, then it is necessarily true that the symmetric matrix 
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dp{dpj 



is positive [40]. In contrast, whenever Eq. (8) is violated, a system with more than one phase 
is energetically favorable. The phase coexistence is governed by the Gibbs' conditions 

<(r,/>r) = /if(r,^f) = ... = <(r,/>r) , ^ = l,...,n, (i4) 

f{T, p-,) = p^iT, p^) = ---= p-iT, pt) , (15) 

where the temperature is the same in all phases. The local positivity conditions on J-'ij give 
rise to a set of n inequalities that divide the parameter space {T, pi} into stable and unstable 
regions. It is important to realize, however, that under conditions of phase separation, Eq. (8) 
may still be valid locally (that is, for ~ Pi)^ but it may nevertheless be possible to 

find significantly different densities p'- and p'- that violate this condition. This leads to the 
existence of metastable states. 

These ideas are illustrated for a simple one-component system in Eig. 1. Here the free 
energy density at fixed temperature is shown as a function of the density. At point D, we 
evidently have 

d'T „ 
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FIG. 1. Free-energy density as a function of density for a one-component system at fixed 
temperature. For clarity, we show T — e^p = T -\- (15.75 MeV)/>, since the variations in T are quite 
small. The additional term linear in the density has no effect on the stability conditions or on 
Gibbs' criteria. 

and the system is stable. At points A and i?, is still locally convex, but these points 
share a common tangent (chemical potential). Thus, although the stability criteria remain 
valid locally around A and i?, the existence of the common tangent implies that these two 
phases with different densities can coexist. The values of T and p at the points A and B 
lie on the binodal surface in the {T, /?} plane, which is a one-dimensional surface in this 
simple example. At point C, the free energy is concave, indicating that this configuration is 
unstable against phase separation. For points between A and B where the second derivative 
of remains positive, the system is metastable and can exist temporarily in a single phase, 
allowing for supercooling or superheating. However, there is an infiection point between A 
and C, and also one between B and C, where the second derivative of vanishes; these 
points lie on the spinodal surface, which delimits the onset of instability. 
According to Gibbs' phase rule [41], at most 

A max — n -\- 2, 

phases can coexist in equilibrium in a system with n conserved charges. Each single phase 
is characterized by the set (T, V") of n + 2 variables, so that if there are K phases, the 
total number of variables is K[n + 1) + 1. (All phases have the same temperature.) We 
assume that the n total charges Qi and the total volume V are also specified, which brings 
the total number of variables to A^(n + 1) + n + 2. These variables are restricted by the Gibbs' 
conditions (14) and (15), which provide [K — l)(n + 1) constraints, and by the conservation 
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laws (11) on the charges (n constraints) and (13) on the total volume (1 constraint). We 
therefore end up with a set of n + 2 independent variables, which we can take to be (T, Qi, V) 
or (T^pi^V). This result implies that the number of degrees of freedom in a system with 
strictly conserved charges is independent of the number of phases [42,43]. 

It is important to emphasize, however, an important difference between a multicompo- 
nent system and one containing a single conserved charge. Although the ratios of the total 
charges Qi/Qj remain fixed once the system has been prepared, the ratios can be different 
in different phases. Moreover, since the independent variables T, Qi, and V determine the 
energetically stable state of the system, it is impossible to impose additional constraints. 
For example, as shown in Ref. [20], one cannot demand that the pressure remain constant 
during a phase transition at fixed temperature, for this would violate either the equilibrium 
conditions or the conservation laws (12). Consequently, the common (vapor) pressures and 
chemical potentials vary when the proportions A" of the phases change. These results will 
be illustrated explicitly for the case of asymmetric nuclear matter in Sees. IV and V. 

A necessary condition for the coexistence of more than two phases is that any pair of 
these phases must be in equilibrium. In addition, as mentioned earlier, a system that is 
stable against separation into two phases is also stable against separation into multiple 
phases. Let us therefore focus on two such phases for a moment, which we will denote with 
a prime and a double-prime. The two sets of densities {/o',/?''} that satisfy Eqs. (14) and 
(15) form a surface in the parameter space {T^pi}; this is the phase separation boundary, 
or binodal. For n conserved charges and two coexisting phases, Gibbs' phase rule implies 
that the binodal is an n-dimensional surface. We can also show that this surface encloses 
all points satisfying 

p, = (1 - \)p', + \p': , < A < 1 , (16) 

that lead to a single (unstable) configuration with a higher value for the free energy. Al- 
though this statement is reasonable from a physical point of view, it is not immediately 
obvious. 

To prove it, we define a tangent plane T^p^j[pi) that is a function of the pi and that can 
be attached to any point in the parameter space pi. Since the slopes at any such point are 
given by the chemical potentials Pi{pi) [see Fq. (7) and Fig. 1], these tangent planes can be 
expressed as 

n n 
T{pO(/?0 = ^ {P^ - Pi)lJn{pi) + ^(r. Pi) = ^ p^P^{pi) " p{T , pi) . (17) 

8 = 8 = 

From the first equality, it follows immediately that when jF(r, pi) is a convex function of the 
Pi, the tangent plane will lie below (the stable region), and when the convexity of is lost 
(that is, when acquires a "saddle point"), the tangent plane will lie above (the unstable 
region) [40]. Moreover, from the second equality, the existence of a common tangent plane 
connecting two distant configurations implies equal pressures and chemical potentials in 
these configurations, which are simply Gibbs' criteria for phase equilibrium. These common 
tangent planes therefore define a binodal surface and allow for phase separation; it also 
follows that every point in the parameter space that leads to an unstable single-phase system 
must lie inside the region enclosed by the binodal surface. This is evident in Fig. 1 for the 
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unstable systems in the neighborhood of point C, which he inside the binodal surface that 
contains A and B. 



-0.8 -0.7 -0.6 




-0.75 [ 

0.05 0.10 0.15 0.20 
P [fm'l 

FIG. 2. Free-energy density as a function of density and asymmetry ps/p for a system with 
two conserved charges at fixed temperature. We again show T — e^p for clarity. The curve for 
symmetric matter has /O3//O = 0. 

In Fig. 2, we show the corresponding situation for a system with two conserved charges 
at hxed temperature. Here p denotes the density of the sum of the charges, and p^ denotes 
the density of the difference; thus, p^j p is a measure of the asymmetry. The dashed curve 
indicates the intersection of the two-dimensional binodal surface obtained from the set of 
common tangent planes with the plane dehned by T = 10 MeV. (The endpoints of the dashed 
curve in Fig. 2 correspond to points A and B in Fig. 1.) Observe that all conhgurations 
where the free-energy density has a saddle point are contained within the binodal. Note also 
that in this example, the free-energy density is always convex with respect to variations in 
/J3 at hxed p and T . 

One feature of the binodal surface is that it may contain critical points. At the criti- 
cal points, if they exist, the two phases can no longer be distinguished by their densities. 
Therefore the critical points form a line that divides the binodal surface into different regions 
describing either a high density (liquid) or a low density (gas) phase. Finally, we note that 
more than two phases can coexist if and only if each pair of phases form a binodal, and if 
all these binodals have a common region of intersection [41]. 

The binodal surface determines the stability boundaries of the system, but it remains 
to show how the system behaves inside, i.e., how to interpolate within the metastable and 
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unstable regions using a Maxwell construction. To explain this in more detail, we consider 
an isothermal compression in a situation where the system can separate into two phases. 



Pa=Pb 



1 = const. 


Binodal / 
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\ / 
\ / 
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1 P3/P = 





Pa Pb 

P 

FIG. 3. The Maxwell construction in a one-component system. The binodal line is obtained 
from similar curves at other temperatures. 

We begin with the familiar case of a one-component system, as shown in Fig. 3. Suppose 
that during the compression, the system encounters the binodal at some point A in {T, p\ 
space. At this point, the whole volume is occupied by a phase with density /j^, and a second 
phase with density is about to emerge in an infinitesimally small volume. The two phases 
at A and B are connected by the Gibbs' conditions, so that they have equal temperatures, 
pressures, and chemical potentials. In a one-component system, B is the point at which 
the system leaves the two-phase region; on the (P, p) diagram of Fig. 3, A and B are 
connected by a horizontal line, the well-known Maxwell construction. In a multicomponent 
system, however, as depicted in Fig. 4, the ratios of the charges in the emerging phase at 
B are generally different from those in the original preparation at A, thus violating the 
conservation laws. The system must therefore evolve instead through configurations that 
maintain the ratios of the total charges (the curve AD), and it leaves the instability region 
at the point D, which lies together with A on the line of constant ratios QijQj. At this 
point, the original phase is present in infinitesimal quantities with densities {/of }, while the 
newly created phase has evolved to point D. The configuration at D is consistent with the 
conservation laws, and in general, the pressure and chemical potentials in the coexisting 
phases have changed throughout the transition. 
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Pa Pc Pd Pb 

P 

FIG. 4. The Maxwell construction in a two-component system. The asymmetry is held constant 
at /O3//O = —0.4 throughout the phase separation. Note that only the solid and dotted curves have 
PsIp = -0.4. 

To determine the nature of the system between these extreme values, we must solve 

p, = (1 - A)^:- + A^:' (18) 

for given values of pi^ with p'- and p'-' lying on the binodal surface. It is important to realize 
that Eqs. (18) are indeed a set of n equations in n unknowns, since the pi are specified, 
and among the 2n + 1 variables and A, n + 1 can be eliminated by virtue of the 

Gibbs' conditions (14) and (15). Moreover, these equations yield solutions with qualitatively 
different characteristics. If the solutions yield all values of A in the interval [0, 1], so that 

=pt for A = , 
p' = pF for A = 1 , 

then the system has undergone a phase transition. However, anticipating the subsequent 
discussion, there are also solutions with < A < A^ax < 1. In this case the system becomes 
unstable to phase separation, but undergoes a retrograde condensation: after occupying a 
maximal volume fraction Amax, the new phase begins to disappear, and the system leaves 
the instability region in the original phase. In either situation, Eq. (18) provides the desired 
Maxwell construction that determines the free energy in the transition region according to 

^(r,^0 = (i- Wr, />:•) + A^(r,^:') . 
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Densities related to other extensive quantities can be computed accordingly. 

We close this section by specializing the general formalism to asymmetric nuclear matter, 
a system of interacting neutrons, protons, and mesons. Such a system is characterized by 
two conserved charges: the total number of baryons 

B = N, + N^ = Vp (19) 

and the total charge, or equivalently, the third component of isospin 

N - N 1 

h = '-^^^^'-Vp,. (20) 

Thus we have 

J^(r, p, p3) = -p{T, /i, /is) + PP + p3p3 ■ (21) 

The stabihty condition (8) imphes the following set of inequalities on the convex free energy 
density: 

> , (22) 



2 



dp 

dp^ dpi \dp3dpj 

To make the physical content of these conditions more transparent, it is convenient to 
introduce the proton fraction y defined by 

y ^ = ^ (24) 

and to rewrite the free-energy density as 

J^iT,p,p3) = nT,P,y) . (25) 
It is straightforward to show that the conditions (22) and (23) are equivalent to 



where we have introduced chemical potentials for protons and neutrons defined by 

Pp = p + ps , pn = p - P3 ■ (28) 

The first inequality is the familiar requirement that the isothermal compressibility is positive, 
that is, that the system is mechanically stable. The second condition refiects the special 
character of the binary system. It expresses "diffusive stability", which guarantees that 
energy is required to change the concentration in a stable system while holding the remaining 
variables (pressure and temperature) fixed [42]. 
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III. RELATIVISTIC MEAN-FIELD EQUATION OF STATE 



To describe the nuclear equation ol state, we use a relativistic mean-field model con- 
taining nucleons, neutral scalar and vector fields, and the isovector p meson field [22,23]. 
The neutral meson fields are self-interacting, including terms through fourth order in the 
fields, as proposed by Bodmer and Price [28,30]. (For simplicity, we omit couplings be- 
tween the scalar and vector fields.) The p meson is introduced in a minimal fashion, as 
discussed in the Introduction. This model can reproduce the observed properties of nuclear 
matter, and recent calculations show that it also gives an accurate description of the bulk 
and single-particle properties of nuclei throughout the periodic table [34]. 

The lagrangian density for this model can be written as 

C = ?[7,(^a^ - g.V^ - \gpT-h^) -{M- g^<P)]i; 
- \ F,^F^' + \ mlV.V^ + ^ (gtiV.V^ f 

-iB,..B-'' + im^b,.b- . (29) 

The scalar, isoscalar- vector, and isovector- vector fields are denoted by ^, V^, and b^, 
respectively, and the vector meson field strengths are F^^ = d^V^ — d^V^ and B^'' = 
d^h"" — d^^h^ — gph^ x b''. We work in natural units with % = c = A;Boitzmann = 1- 

In the mean-field approximation, the pressure p and the energy density £ are easy to 
compute and can be written as [22,23] 

P = ^ [i?5(.„ M*) + i7,(.., M*)] + ^ + ^ + ^ 

L$2 ^3 ^4 /ggX 

2gr e^r 245r ' ^ ^ 

£ = ^ [Hs{p,, M*) + H,{v^, M*) + M*'Hs{pp, M*) + M*'Hs{Pn, M*)] 



+ Wp-^W'-^WU- g.Rps -^R' 
' 2gl 24 2 "' 2gl 



,2 ^ 



777 , , 

2gl 6^3 24^4 



with the conserved baryon density 



p = \ \g^{v,, M*) + Gs{Pn, M*)] (32) 



7r2 



and isospin density 



P^ = j^ [Gsi^p^ M*) - Gs{Pn, M*)] . (33) 
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Here, following Bodmer [30], we define the scaled meson fields $ = g^cf), W = (/yVo, and 
R = (/p^o, with 60 the timelike, neutral part of the p meson field. 

The baryon effective mass and effective chemical potentials are defined in terms of the 
meson mean fields as 

M* = M - $ , (34) 
//, = /i, - 11/ - i i? , (35) 

u^ = p^-W + . (36) 

We also define the required integrals over the thermal distribution functions as 

1 / 1 1 \ 



Vl + e'^I^e^'^)-^] 1 _^ e'^[^('='^)+^] 

k'^-^dk (I I \ 

where E[k^M) = [k'^ + Af^)^/^, and n > to ensure convergence is understood. Note the 
important sign differences in these two relations. These functions obey the useful recursion 
relations (valid for n > 1) 



dG 



- nH.,+, + {n- l)M^H^-i , (39) 
-(n-l)MG'„_i , (40) 



dGn+i 
dM 

dHn+i 
dji 

dHn+i 
dM 



{n - l)Gn-i , (41) 
-{n - l)Mi7„_i . (42) 



Thermodynamic equilibrium requires that the thermodynamic potential be stationary 
with respect to changes in the mean fields, which leads to the self-consistency equations 

H'- f 1 + 4 ^ =4''. (11) 

where the scalar density is given by 

M* 



R = 7r^P3. (45) 



Ps 



[H:,{v,,M*) + H:,{v^,M*)\ . (46) 



These equations allow the fields to be held fixed when computing thermodynamic quantities 
as derivatives of the thermodynamic potential. Moreover, they allow the energy density to 
be expressed as 
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^ = ^ [H^i^p^ M*) + H,{v^, M*) + M*'Hs{pp, M*) + M*'Hs{Pn, M*)] 

+ ^^^ + |^^ + -^/>^ + f|$2 + -^$3 + ^$4, (47) 
2g^ 8 Smj 2g^ 6g-^ 2ig^ 

where $ and W are understood to satisfy the equations given above. 

Finally, we exhibit the corresponding expression for the entropy density, which follows 
from the Gibbs' relation: 

a = {p + £ - fipPp - finPn)/T 

= [mi'^-p. M*) + m,{v^, M*) 

+ 3M*\Hs{pp, M*) + Hs{Pn, M*) 
-3//,G'3(/^„M*)-3//„G'3(/^n,M*)| . (48) 



Note that none of the previous relations involve thermal contributions from the mesons, 
since their masses are too large for these to be relevant. Indeed, at the temperatures of 
interest in this work, the antibaryon contributions are negligible as well. 

To specify the parameters, we observe that models that successfully reproduce bulk and 
single-particle properties of finite nuclei share characteristic properties in infinite nuclear 
matter [32,44]. After taking the zero-temperature limit of the preceding results, one can 
obtain an explicit set of transcendental equations that determines the parameters for the 
desired choice of nuclear matter properties. (We will not exhibit these equations here; see 
Refs. [30,34].) The parameters so obtained are listed in Table I. 



TABLE I. Mean- Field Parameters 
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374.77 


260.57 


106.91 


3.0809 


8.106 


0.02364 



TABLE IF Nuclear Matter Properties 



Equilibrium Properties: 










kp p 


M*/M 
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a4 



1.30fm-i 0.1484fm-3 0.60 -15.75 MeV 250 MeV 35MeV 



Critical Values: 






Tc Pc 


M*/M 


Pc 


14.40 MeV 0.04661fm-3 


0.8543 


0.2010 MeV/fm^ 
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Note that the nucleoli and vector meson masses are chosen to take their empirical values 
[M = 939 MeV, niy = m^^ = 783 MeV, nip = 770 MeV), and only the ratios of couplings to 
masses (denoted by Cf = gf M'^ / mf) are needed in infinite matter. The resulting properties 
of nuclear matter, as well as the properties at the critical point in symmetric matter, are given 
in Table II. To generate acceptable bulk nuclear properties, it is important to accurately 
reproduce the nuclear matter equilibrium density, energy /nucleon eo, baryon effective mass 
Af*, compressibility l/Ky [45], and bulk symmetry energy 04. The first three of these are 
tightly constrained [32], whereas the latter two are not. We will begin by studying warm 
nuclear matter for the values of Ky and 04 given in Table II and later examine the sensitivity 
to reasonable variations in these values. The observant reader will notice that the model has 
six free parameters (in nuclear matter) that are determined by only five constraints. Thus 
there are actually an infinite number of parameter sets that will reproduce the equilibrium 
properties listed in Table II. These sets differ in the way the nonlinear meson interactions 
are split between the scalar terms (/€, A) and the vector term [(). To arrive at the parameter 
values in Table I, Dirac-Hartree calculations of finite nuclei were also performed in this model 
[34], and the parameters were tuned to give optimal bulk- and surface-energy systematics.^ 
This allows us to proceed with the most realistic mean-field nuclear equation of state possible. 

We close this section with some remarks concerning the numerical procedures. Although 
the solution of the self-consistent equations (43)-(45) is in principle straightforward, the 
analysis becomes involved due to multiple roots. For example, at low temperatures, a given 
set {z/p, Pn} leads to either one or three solutions for the scalar field in Eq. (43). We resolve 
this problem using the crucial fact that any quantity can be uniquely and continuously 
parametrized in terms of the effective mass M* = Af — $. To give a concrete example, 
consider the equation of state at constant pressure and proton fraction. In this case we need 
to solve three equations, namely, Eq. (30) with a given value of p on the left-hand side, an 
equation that fixes the proton concentration j/, 

^ Gs{pp,M*) + Gs{Pn,M*) ' 

and Eq. (43). This set is solved for a given value of the effective mass leading to a unique 
root of the form {z/p(Af*), z/„(Af*), r(Af*)}, which can in turn be used to evaluate all the 
remaining quantities of interest, e.^., the density and entropy. Once the solutions have been 
obtained for a given Af*, we can proceed to map out all the desired variables by making 
small incremental changes in this parameter. 

The main ingredient in our thermodynamic treatment is the binodal surface, namely, 
the collection of points in parameter space that satisfy the Gibbs' conditions (14) and (15). 
Numerically this surface is most easily parametrized in terms of the pressure, so that the 
equations we solve simultaneously are 

p = p(//;,//;,M*') , (49) 

p = K<,<',M*") , (50) 



^The authors are grateful to R. J. Furnstahl for his assistance in obtaining the parameter sets. 
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for given values of p and T. This procedure works well except near the critical points, where 
the two solutions coalesce. Correspondingly, results in this region must be obtained by 
interpolation, after one determines the location of the critical points using Eq. (54), below. 
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IV. APPLICATION TO ASYMMETRIC NUCLEAR MATTER 



Phase transitions in binary systems are more complex than in one-component systems. 
In the case ol nuclear matter, as discussed in Sec. II, the global instability boundary forms 
a two-dimensional surface in (T, p, j/) space, enclosing the region where either mechanical 
instability [Eq. (26)] or diffusive instability [Eq. (27)] occurs. 
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FIG. 5. Pressure as a function of baryon density at fixed temperature for various proton frac- 
tions y. 

To be more specific, Eig. 5 shows the pressure as function of the baryon density at a fixed 
temperature but for different proton fractions, using the parameter set in Table I. Eor small 
J/, and in particular, for pure neutron matter (j/ = 0), the pressure increases monotonically, 
so the matter is stable at all densities. In contrast, for y ^ 0.2, the compressibility becomes 
negative, indicating a mechanical instability. The full complexity of the binary system is 
indicated in Eig. 6, where chemical potential isobars for neutrons and protons are shown as a 
function of y at fixed temperature. Above a certain critical pressure the matter is stable, 
but for p < the second condition (27) is violated, and the system becomes chemically 
unstable. The critical isobar pc is determined by the existence of an infiection point: 

This isobar marks the upper boundary of instability with respect to the pressure and defines 
a critical point (pc,J/c) for a given temperature [42]. 
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FIG. 6. Chemical potential isobars at fixed temperature as a function of y. The curves labeled 
a through e have pressures p = 0.25, 0.198, 0.15, 0.10, 0.075 MeV/fm"^, respectively. The curves 
labeled b are at the critical pressure pc = 0.198 MeV/fm"^. 

The Gibbs' conditions (14) and (15) for phase equilibrium demand equal pressure and 
chemical potentials for two phases with different concentrations. Thus the two desired 
solutions form the edges of a rectangle and can be found by means of the geometrical 
construction shown in Fig. 7 [19]. The collection of all such pairs yi[T,p) and y2{T^p) form 
the binodal surface. For a given temperature, the two-phase region is limited from below 
by the pressure at equal concentrations peq, which corresponds to symmetric nuclear matter 
(y = 0.5). Correspondingly, the rectangle in Fig. 7 shrinks to a point at this particular 
pressure, and the system becomes stable again for p < p^q. We did not find any cases of 
three-phase coexistence, which would require two rectangles with a common side. 
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T = 10MeV p = 0.1 MeV/W 
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FIG. 7. Geometrical construction used to obtain the proton fractions and chemical potentials 
in the two coexisting phases at fixed T and p. 




FIG. 8. The binodal surface indicating the two-dimensional phase-coexistence boundary is 
shown in (p,T,y) space. The critical temperature Tc(y = 0.5), the line of equal concentrations 
(LEG), the line of critical points (LCP), and the line of maximal asymmetry (LMA) are indicated. 
T is in MeV, p is in MeV/fm"^, and y is dimensionless. 
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The phase-separation boundary, or binodal surface, obtained from the preceding geo- 
metrical constructions [see Eqs. (49) through (53)] is indicated in Fig. 8. (The shape of this 
surface has been previously described as a "filet mignon" [19].) Several slices at constant 
T are indicated, and one observes that the enclosed area in these sections decreases with 
increasing temperature until it vanishes at the critical point Tc of symmetric nuclear matter. 
This is also the point at which the lines LCP and LMA meet at y = 0.5. The line of critical 
points (LCP) is determined by the solutions of Eq. (54), and it can be parametrized uniquely 
by just one of the coordinates (r,p, j/). The LCP begins at the critical point of symmetric 
nuclear matter [y = 0.5, T = Tc) and ends at j/ = j/o = 0.057 at T = 0. Note also that for 
a given temperature, the critical point determines the maximum pressure in the two phase 
region. 

In addition to the LCP, we have also indicated the points on the binodal surface with the 
maximal asymmetry (LMA), or minimal proton fraction j/min, at each temperature.^ Any 
system with y < ymin is external to the two-phase region and is therefore stable. Moreover, 
since ymin is a monotonically increasing function of the temperature, the LMA determines 
the maximum temperature T^ax of phase separation at any given y [19]. The LMA also 
begins at the critical point of symmetric nuclear matter and descends to a small, positive 
value of j/inin at r = 0. This implies that pure neutron matter is stable at all temperatures 
in this model. 

Configurations that separate into two phases each having equal numbers of neutrons and 
protons form the line of equal concentration (LEC). In these cases, the binodal section at 
any T degenerates into a point, at which the relation [41,43] 



is satisfied. (This is an example of the Gibbs-Konowalow rule [43].) The LEC coincides 
with the projection y = 0.5 in our model, which implies that only symmetric nuclear matter 
will separate into phases with equal concentration [19]. Since symmetric matter behaves 
as a one-component system,"^ one observes how the phase separation simplifies in this case: 
in an isothermal compression, the system evolves until it encounters the binodal (which 
is just the LEC) and then remains there until the transition is complete. In contrast, for 
y ^ 0.5, the system encounters a two-dimensional section of the binodal surface. Since the 
energy contains a term proportional to [see Eq. (47)], it is energetically favorable for 
asymmetric matter to separate into a liquid phase that is less asymmetric and a gas phase 



^As with all results in this section, our analysis is symmetric with respect to protons and neutrons, 
so that only the physically useful regime < y < 0.5 need be considered. 

■^Strictly speaking, symmetric nuclear matter is an azeotrope. During the liquid-gas phase sepa- 
ration, all equilibrium configurations have two phases with the same composition (y = 0.5). This 
type of phase equilibrium is called "indifferent equilibrium" [43]. Since the liquid-gas phase equi- 
librium in a one-component system is also indifferent, symmetric nuclear matter behaves as a 
one-component system. 




(55) 
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that is more asymmetric, rather than into two phases with equal concentration. This leads 
to more complex phase separations, as we discuss shortly. 

For a given temperature, the binodal section is divided into two branches by the critical 
point and the point of equal concentration. One branch describes the system in a high- 
density (liquid) phase, while the other branch describes the low-density (gas) phase. These 
two branches contain the beginning and ending configurations of the phase transition. 

We now return to the behavior of the matter under isothermal compression, to illustrate 
the different phase-separation scenarios. 
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FIG. 9. Binodal section at T = 10 MeV. The points A through D denote phases participating 
in a normal phase transition. The critical point (CP) and the points of equal concentration (EC) 
and maximal asymmetry (MA) are also indicated. 

Consider the situation in Fig. 9, which shows a section through the binodal surface at 
T = 10 MeV. The critical point CP, the point of maximal asymmetry MA, and the point 
of equal concentration FC are indicated, and the validity of Fq. (55) is apparent. Assume 
that the system is initially prepared in the low-density (gas) phase with proton fraction 
y = 0.3. During the compression, the two-phase region is encountered at the point A, and 
now a (liquid) phase with a higher density begins to emerge. The geometrical construction 
described above determines the density and the proton fraction yg of this new phase, which 
occurs at the point labeled B. As the system is compressed, the total proton fraction y 
remains fixed, as dictated by the conservation laws, but the gas phase evolves from A to D, 
while the liquid phase evolves from B to C . At the point C, the system leaves the region of 
instability. The original (gas) phase is about to disappear, and it exists in an infinitesimal 
volume with a density and proton fraction y^ corresponding to the point D. 
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FIG. 10. Isotherms for a normal phase transition at T = 10 MeV and initial condition y = 0.3. 
The Maxwell construction produces the curve AC. Note that 0.3 yg. 

To evolve the system between configurations A and C, we must solve the equations 

/> = (!- + A/ , (56) 
p, = {l- X)p', + A^;; , (57) 

according to (18), for densities that lie on the binodal surface, and for given values of p and 
P3 = [2y — l)p. The result is the generalized Maxwell construction in the binary system. 
The corresponding isotherms are drawn in Fig. 10. The dotted line between A and C 
is the unphysical course of the pressure at the fixed total proton fraction, and the nearly 
straight line connecting A and C is the interpolation due to the Maxwell construction, which 
corresponds to the stable configuration at each intermediate density. As mentioned earlier, 
the compression in the two-phase region is nonzero because of the pressure change forced by 
the conservation laws, in contrast to the behavior found in a single-component system. The 
volume fraction A starts with A = at A and runs through the whole interval [0, 1], ending 
with A = 1 at C. Since the points A and C lie on different branches of the binodal surface 
(as defined above), the matter has undergone a phase transition from a gas to a liquid phase. 

Interestingly enough, the geometry of the binodal surface offers a second possibility. In 
the previous example, there is a transition between the two branches of the binodal surface 
because the value of y in the original phase is larger than j/c- For y < yc, however, the 
system enters and leaves the two-phase region on the same branch, so the system remains 
in the same phase. 
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FIG. 11. Binodal section at T = 10 MeV. The points A through D denote phases participating 
in a retrograde phase transition. 

This situation is depicted in Fig. 11. We consider the T = 10 MeV isotherm and prepare 
the system with j/ = 0.15 < j/c- The system becomes unstable at the point A, and as before, 
a hquid phase with a higher density begins to emerge at B. The system is compressed at 
fixed total j/, with the hquid phase evolving from B to D, and the gaseous phase, from A to 
C . At C the system crosses the binodal again, but this time on the same branch, that is, still 
in the original (gas) phase. The high-density (liquid) phase, now at D, vanishes at this stage. 
The Maxwell construction for the corresponding isotherms, which follows from the solution 
of Eqs. (56) and (57), is represented by the solid line connecting A and C in Fig. 12. Note 
that this new phenomenon is caused solely by a diffusive instability. The matter remains 
mechanically stable throughout the entire process, as indicated by the dotted curve, but 
it is energetically favorable to separate into two phases with different proton fractions. In 
contrast to the previous case, one also finds a different behavior for A. The initial value 
A = at A increases up to a maximal value A^ax < 1 and then decreases to zero when C 
is reached. Thus, although a second phase is present between A and C, the system does 
not convert completely; on the contrary, the second phase vanishes steadily after having 
occupied a maximum volume fraction {Xmax)- This retrograde condensation is unique to the 
binary system and does not occur in one-component systems. 
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FIG. 12. Isotherms for a retrograde phase transition at T = 10 MeV and initial condition 
y = 0.15. The Maxwell construction produces the curve AC . 

A more complete picture of the different Maxwell constructions is shown in Fig. 13. In 
part (a), we consider different proton fractions at a fixed temperature. At j/ = 0.5 (symmetric 
matter), we obtain the familiar result with a constant vapor pressure, represented by a 
horizontal line. With increasing asymmetry (decreasing j/), the compression increases in the 
two-phase region, as the vapor pressure is no longer constant. For j/min < j/ < j/c, the matter is 
mechanically stable, and the system undergoes retrograde condensation. Finally, the system 
becomes completely stable for y < j/min- Part (b) shows similar Maxwell constructions on 
different isotherms at fixed proton fraction y = 0.4. For a more thorough discussion of the 
variation of the densities in the two phases throughout the transition, see Ref. [19]. 
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FIG. 13. (a) Maxwell constructions at fixed temperature for various proton fractions, 
(b) Maxwell constructions at fixed proton fraction for various temperatures. 
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To close this section, we consider the order of the hquid-gas phase transition. First- 
order phase transitions in single-component systems are characterized by discontinuities in 
certain physical quantities. This occurs because quantities like the density and the entropy 
are different in the two distinct phases, but they remain constant throughout the phase 
transition. This behavior does not occur in the binary system, because the constraints of 
charge conservation and Gibbs' criteria force the density (and pressure) in each individual 
phase to change throughout the transition. We might therefore expect the transition in the 
binary system to be "smoother". To make this point more precise, we consider the Gibbs 
free energy (or free enthalpy) per nucleon 



G{T,p,N,,N^ 



9iT,p,y) = y^p + (1 - y)^r 



(58) 



which allows us to discuss isobaric processes as a function of temperature. 

This quantity is shown for several different asymmetries in Fig. f4. As discussed in 
Sec. II, the Gibbs free energy must be computed according to 



9{T, p,y) = [1- X)g{T, p, y') + \g{T, p, y") 



(59) 



in the transition region, where A now specifies the number of particles in each phase: 

iV" + N" 

X = ^^—^ , (60) 



and where y' and y" denote the different proton fractions. If one heats the symmetric system 
through the transition point, as indicated in Fig. I4a, it remains at constant temperature 
until the transition is completed, producing a striking kink in the free enthalpy curve. 
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FIG. 14. Specific free enthalpy at constant pressure as a function of temperature, (a) Free 
enthalpy in symmetric matter. The solid curves represent the stable single-phase configurations, 
and their dotted extensions denote metastable systems. The matter is unstable along the remaining 
dotted curve, (b) Same as (a) for y = 0.3, showing that the free enthalpy is continuous throughout 
the phase transition. The points Ta and Tg denote the endpoints of the coexistence region, (c) Free 
enthalpy for various proton fractions. The values are y = 0.5, 0.45, 0.4, 0.35, 0.3 for curves a through 
e, respectively. 
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In contrast, for matter with Np ^ Nn (Fig. 14b), the phase transition begins at a 
temperature Ta and ends at Tb > Ta, leading to a completely smooth free enthalpy curve. 
The free enthalpy for various values of the proton fraction is shown in Fig. 14c to illustrate 
how the kink in g develops as the system becomes more symmetric. 

It is also of interest to examine the entropy per nucleon, which is defined by 



,,(r.p.,) = | = -(M|M) 



(61) 



The entropy in the transition region requires special attention, because only the total proton 
fraction is held constant in Eq. (61). Therefore 

s{T,P,y) = (1 - X)s{T,p,y') + Xs{T,p,y") + (^j [9{T,p,y') - g{T,p,y")] 



(1-A) 



'dl\ f dgjT, 



p,y 



dT / V dy' 

/ p,y \ 



'df_\ f dgjT, 



p,y 



T,p 



dT / \ dy' 

/ p,y \ 



T,p 



(1 - \)s{T,p,y') + \s{T,p,y") 



(62) 



where the last identity can be verified by observing that y = {1 — X)y' + \y" and that the 
chemical potentials are identical in both phases. Note that in evaluating higher derivatives, 
e.g., the heat capacity, similar care must be used to include the temperature dependence of 
A, J/', and y" . 
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FIG. 15. Specific entropy as a function of temperature at constant pressure for various asym- 
metries. The curves have the same values of y as in Fig. 14c, with y = 0.5 at the bottom. 
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Figure 15 shows the specific entropy for several different asymmetries. The kink in g for 
the symmetric system gives rise to a discontinuity. In contrast, the entropy in an asymmetric 
system is continuous, but there are now kinks at the end points of the phase transition. By 
using 



we can calculate the amount of heat that has to be transferred to the system during the 
course of the transition. For the symmetric system, this yields the well-known expression 
for the latent heat: 



which means that the all of the energy is used to convert the liquid into vapor. 

The behavior of asymmetric matter is quite different. The temperature changes during 
the transition and so does the entropy, so that the integration becomes nontrivial: 



and thus some fraction of the energy is used just to heat the system, as is evident from 
Fig. 15. The concept of latent heat is therefore not strictly applicable to the transition in 
the binary system. A well-known example of this behavior is the distillation of alcoholic 
beverages; the concentration of alcohol in the liquid changes throughout the distillation, 
which produces a change in the boiling point. 



dQ = Tds , 



(63) 



Ql = T{sa - Sb) , 
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FIG. 16. Specific entropy as a function of temperature for different isobars. (a) Sym- 
metric nuclear matter with y = 0.5. The pressures on the curves labeled a through e 
are p = 0.05, 0.10, 0.15, 0.20, 0.25 MeV/fm"^, respectively. Curve d is at the critical pressure 
p = = 0.20 MeV/fm"^. (b) Asymmetric nuclear matter with y = 0.3. The curves are labeled 
as in (a), and the critical pressure (curve d) is still pc = 0.20MeV/fm"^ (to two significant figures). 

Figures 16a and 16b show the entropy as a function of temperature for several different 
isobars in symmetric and asymmetric matter, respectively. By examining the curves at the 
critical pressure one observes that [ds/ dT)p^y becomes infinite in symmetric matter, while 
it remains finite and positive in the asymmetric case. 

Finally, the distinct behavior of the binary system has an even more dramatic impact on 
the heat capacity per nucleon: 




which is illustrated in Fig. 17. In symmetric matter, the discontinuity in the entropy pro- 
duces an undefined heat capacity at the transition temperature. In contrast, in the general 
case, there are finite discontinuities in the heat capacity at the endpoints of the transition 
region, with finite values of Cp in between. Thus, according to Ehrenfest's definition of phase 
transitions [46], the first- order phase transition in symmetric matter becomes a second-order 
transition in the asymmetric case. In particular, in the present model, the phase transition 
occurs over ~ 1.5 MeV for matter that is 40% protons and over ~ 5 MeV for matter that 
is 30% protons. The latter concentration might be obtained in energetic collisions with 
radioactive ion beams. 
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FIG. 17. Specific heat capacity as a function of temperature for several proton fractions. For 
clarity, the curves for asymmetric matter have been scaled by the factors indicated in parentheses. 
Note the logarithmic scale on the left. 
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V. HEAVY-ION COLLISIONS 



We consider the energetic collision of two heavy ions within our simple thermodynamic, 
hydrodynamic, and mean-field picture, and concentrate on the new features that arise as a 
function of proton fraction. Although our discussion of the theory will encompass the entire 
regime of y < 0.5, to obtain reasonable estimates for the empirical size of the new effects, 
we will restrict consideration to 0.3 < j/ < 0.5, where the lower value might be obtainable 
with radioactive ion beams. We shall assume that the combined system is compressed and 
heated and ultimately reaches equilibrium at some finite temperature, density, and pressure. 
As is well known, the question of whether the system actually reaches such an equilibrium 
state is a difficult one, which we will not attempt to answer here. We simply assume that 
such a state arises and follow the subsequent expansion and cooling of the nuclear matter. 
The properties of the system can be deduced from Figs, f 8 through 20, where the pressure 
is shown as a function of the baryon density for various proton fractions, temperatures, 
and specific entropies. Also indicated are the coexistence curves or binodals, which are 
determined by the Maxwell construction discussed earlier. 
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FIG. 18. Pressure as a function of density at fixed temperature for various y. The critical point 
(CP), coexistence curve (CE), diffusive spinodal (DS), and isothermal spinodal (ITS) are indicated. 
The dashed curves are discussed in the text. 
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FIG. 19. Pressure as a function of density at fixed y = 0.3. The solid curves are labeled by the 
temperature (in MeV). The critical point (CP), coexistence curve (CE), diffusive spinodal (DS), 
and isothermal spinodal (ITS) are indicated. The dashed curves are discussed in the text. 

The spinodals that determine the boundaries of the unstable region now generally come 
in three varieties, one arising from a mechanical instability, one arising from a thermal 
instability, and one arising from a diffusive (or chemical) instability. The conditions that 
determine the different spinodals will be discussed shortly. The region between the binodal 
and the most extensive spinodal contains metastable states, which correspond either to 
superheated liquid or supercooled (supersaturated) vapor. 

The metastable states play a key role in understanding the phase transition. In stable 
configurations, Eq. (8) is satisfied, so that 



In the metastable region, two phases can be found that produce AJ-;, < 0, but the energy 
barrier A J-;, remains positive for small variations in p and y about the single-phase values. 
Thus the free-energy density is locally convex, and the energy barrier prevents phase sep- 
aration due to infinitesimal fiuctuations. Eventually, A J-;, becomes negative as the system 
enters the labile region, signaling the instability even to infinitesimal fiuctuations. We post- 
pone a discussion of the various mechanisms that can produce the phase separation and 
concentrate here on the behavior of the free energy in the region of the spinodals. 




(65) 



p = {l- \)p' + \p 




(66) 
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FIG. 20. Pressure as a function of density at fixed y = 0.3. The solid curves are labeled by 
the specific entropy. The critical point (CP), coexistence curve (CE), diffusive spinodal (DS), 
isothermal spinodal (ITS), and adiabatic spinodal (AS) are indicated. The dashed curves are 
discussed in the text. 



According to the discussion in Sec. II, we have 

( dp\ 



AJ^b > for 



> and 



T,y 



> 



(67) 



in the stable and metastable region. The content of these relations becomes more transparent 
if we consider small fluctuations around the equilibrium density and concentration: 



P' 



p' = P + ^P 

II I A 

- P + Ap 



II II 

, y 



y + ^y' ■, 
y + Ay" . 



(68) 



An expansion through second order in small quantities produces 

' dpp\ ( dpr. 



2A 



dy 



+2p 



'dpp 
. dp , 



T,y 



. dp 



dy 

Ay'Ap'] + ... 



I A „/ 



'Ay 



l\2 



(69) 



where Ap" {Ay") has been eliminated in favor of A and Ap' {Ay') using Eq. (66). Consistent 
with our general discussion, this bilinear form will be positive if 
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(70) 



and if 




, dy 




(71) 



Therefore, as long as the conditions (70) and (71) are valid, the system remains on a parabolic 
free-energy surface, with a minimum at Ap = Ay = 0. The spinodals are defined by the 
densities at which these inequalities become invalid. As is clear from Fig. 18, there are 
generally four such densities of interest in an isothermal expansion, corresponding to the 
intersections of the isotherm with the spinodals. Each of these four densities marks a change 
in the sign of the derivative [dfip/ dy)T,p- At the largest density, this derivative changes sign 
from positive to negative, indicating the onset of the diffusive instability; nevertheless, the 
compressibility remains positive, so that only the second inequality (71) is violated. This 
also implies that the free-energy density now has a saddle point, being stable against density 
fiuctuations at fixed j/, but unstable against fiuctuations in concentration at fixed pressure. 
Through continued expansion, one next encounters the mechanical instability, at which 
point [dp/ dp)T,y becomes negative, but [dpp/ dy)T,p now becomes positive again, so that 
both conditions (70) and (71) are violated. Thus the free-energy surface still has a saddle 
point, but it has rotated in the (/?, y) plane. As the density decreases further, one eventually 
reaches densities where the derivatives change sign in reverse order, producing the remaining 
two intersections of the spinodals with the isotherm. 

The conclusion from this analysis is that the spinodal structure is qualitatively different 
in the asymmetric system than in the symmetric one. There are now two types of fiuctu- 
ations: one corresponding to changes in p (isoscalar) and one corresponding to changes in 
y (isovector). Due to the form of the nuclear symmetry energy, the diffusive spinodal (DS) 
encloses more of the configuration space than does the isothermal spinodal (ITS). The DS 
passes through the critical point and includes all densities where there is also a mechanical 
instability, which occur in the region bounded by the ITS. Therefore the diffusive instability 
defines the relevant spinodal for the asymmetric system. The mechanical instabilities are 
restricted to isotherms with y > y^^ where j/* is defined by the infiection point 




T,y* 



dp\ 



. 



T,y* 



As expected, both spinodals coincide at y = 0.5, where only mechanical instability is possi- 
ble. Note also that there are isotherms that pass through the metastable region and never 
intersect either spinodal. These isotherms, along which retrograde condensation is possible, 
never become labile and allow for the system to evolve completely through the metastable 
region, if the process is carried out carefully enough (or fast enough). 

So far, we have considered only isothermal processes, but in fact, the equilibrium evolu- 
tion of the system will depend strongly on the variables that actually remain fixed during 
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the process. As we have noted, it is not clear that any thermodynamic variables can be 
assigned to realistic situations encountered in warm, expanding nuclei. Nevertheless, in- 
tranuclear cascade calculations [8,14] suggest that this expansion will be isentropic. To 
discuss isentropic (adiabatic) processes, we must generalize our formalism slightly. The ap- 
propriate state function to describe trajectories at constant specific entropy s = S/ B is the 
energy/baryon, E/B = e(s,u,j/), where v = l/p. The corresponding stability condition can 
be formulated as: 

e(5, V, y) <{1- A)e(s', v' , y') + Ae(s", v'\ y") , (72) 

with 

s = {l-X)s' + \s\ ^; = (l-A)^;' + A^;", y = (1 - A)j/' + V , (73) 

where A [see Eq. (60)] determines the fraction of particles in each phase. 

As in the case of the free-energy density, the global criterion Eq. (72) can be reformulated 
in terms of local conditions, namely, 

(74) 
(75) 
(76) 



1 / dp^ 
c.p \dp^ 



s,y 




where Cp is defined in Eq. (64). The first two inequalities are relevant for both symmetric 
(one-component) and asymmetric (two-component) systems, and the third embodies the 
diffusive or chemical stability criterion that arises for asymmetric matter. Thus there will 
generally be three spinodals, as indicated in Eig. 20. 

To connect this to our earlier discussion of isothermal processes, we use the relation [47] 




s,y 




(77) 



to rewrite Eqs. (75) and (76) as 




(78) 
(79) 

The equivalence of these two ways of writing the stability criteria implies that the stability 
boundaries are independent of the actual process. Moreover, since in our model (see also 
Ref. [19]), Cy = T[ds/ dT)y^y is always positive, the conditions (78) and (79) reduce to those 
discussed earlier in Eqs. (70) and (71). Thus, phase separations in isothermal processes 
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in symmetric matter occur only due to mechanical instabilities (density fluctuations), as 
determined by the ITS. In asymmetric systems, the final stability condition is also relevant, 
leading to diffusive instabilities and the appearance of the DS as well. 

In contrast, Cp undergoes several sign changes in an adiabatic expansion, so that all three 
criteria (74)-(76) must be considered. One sign change accompanies a sign change in the 
isothermal compressibility [or [dp/ dp)T,y]^ and one accompanies a corresponding change in 
the adiabatic compressibility [or [dp/ dp)s,y]- The resulting situation is depicted in Fig. 20. 
The outermost spinodal (DS) defines where Eq. (76) is violated, which occurs due to the 
diffusive instability, since both c,p and the compressibilities are all positive here. The middle 
spinodal is determined by condition (75), which is violated when c,p becomes negative; be- 
cause of Eq. (77), this curve is identical to the ITS, since Cy and [dp/ dp)s,y are still positive, 
so [dp/ dp)T,y must also become negative. In other words, the ITS encloses the region of 
mechanical instability in an isothermal expansion and also encloses the region of thermal in- 
stability in an adiabatic expansion. The mechanism for this thermal instability is interesting, 
because Cp changes sign by passing through infinity rather than zero. Thus a small thermal 
fiuctuation at constant pressure near the ITS creates a density fiuctuation, but leaves the 
temperature unchanged (since Cp is very large); since the isothermal compressibility becomes 
negative, the system is unstable to the induced density fiuctuation. 

Einally, the innermost curve is the adiabatic spinodal (AS), which determines when 
Eq. (74) is violated, signaling the onset of mechanical instability in the adiabatic process. 
Here Cp becomes positive again by Eq. (77), since both compressibilities are negative here. 
Note also that, as discussed earlier, [dpp/ dy)T,p changes sign from positive to negative at the 
DS and then changes back again at the ITS, so that once any stability criterion (74)-(76) 
is violated, it remains invalid as one proceeds deeper into the unstable region. 

Just as we found in the case of isothermal processes, the DS separates the metastable 
and labile regions in isentropic processes. Configurations that indicate sign changes of the 
compressibilities and of Cp are included by this boundary. We also observe that even in 
symmetric nuclear matter, the ITS determines the boundary of the labile region in both 
isothermal and isentropic processes, which has often been overlooked in previous studies 
[7,8,11,12,36,48]. Although this result apparently implies that the AS is irrelevant, this 
conclusion is premature, since the actual mechanism for spinodal decomposition depends 
on the relative rate of thermal and mechanical fiuctuations as the warm matter expands 
through the spinodal region. 

Einally, because of the numerous relations between the thermodynamic variables and 
their derivatives enumerated above, correct results will be obtained only if the underlying 
calculation of the nuclear equation of state is thermodynamically consistent. This is known 
to be true for the relativistic MET used here [22,36], but consistency is difficult to maintain 
in more sophisticated approximations to the relativistic many-body problem [49]. 

Armed with this understanding of the spinodal structure, we turn now to the evolution 
of the warm matter. 



"^Note that the usual definition of the isothermal compressibility [47], = —[dV / dp)T/V can be 
rewritten as p[dp/ dp)T = 1/kt- Analogous expressions hold at fixed entropy. 
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FIG. 21. Properties of nuclear matter as functions of temperature and density. The dotted 
curves are contours of equal energy/baryon (in MeV), and adiabats are shown as dot-dashed curves. 
The solid lines denote the coexistence curve (CE), the diffusive spinodal (DS), the isothermal 
spinodal (ITS), and the adiabatic spinodal (AS). Part (a) shows results for asymmetric matter 
with y = 0.3 and part (b), for symmetric matter {y = 0.5). 
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The relation between the evolution and the thermodynamic variables can be studied 
using Fig. 21, which shows the properties of the matter for this equation of state in the T, p 
plane for two different values of y. Equipotential surfaces and adiabats are also indicated. 
The relations used to determine the nuclear matter properties at very low densities are given 
in the Appendix. 

The solid curves in the lower-left corners of Figs. 21a, b determine the phase behavior. 
The outer curve labeled CF is the coexistence curve, which is simply a section through the 
binodal surface at fixed y = 0.3, 0.5. Inside this curve, the stable configuration is a mixture 
of liquid and gas. Also shown is the diffusive spinodal (DS), which is determined by the 
condition [dfip/ dy)T,p = 0, as discussed earlier. Between the spinodal and the CF, the 
system can exist in metastable superheated or supersaturated states. 

There are several qualitative differences between the curves for asymmetric matter and 
for symmetric matter (see also Fig. 8 in Ref. [36]). First, the most extensive spinodal is 
defined by the diffusive instability rather than by thermal or mechanical instabilities. (We 
note that mechanical instability is often assumed to be the relevant one even in asymmetric 
systems.) Second, as is evident from Fig. 8, the critical temperature Tc{y) is not unique, 
but varies with the concentration along the LCP. Finally, the maximum temperature of 
phase separation does not occur at the critical point, but rather along the LMA; it is 
apparent from Fig. 21a that this occurs at a lower density than the density at the critical 
point. We observe, however, that although these differences exist in principle, for physically 
accessible systems, the temperature differences are small. For example, in the present model, 
the critical temperature changes from rc(0.5) = 14.4 MeV to rc(0.3) = 13.1 MeV, and 
rmax(0.3) — rc(0.3) ~ 0.2 MeV. Since statistical fiuctuations near the critical point are 
expected to prevent the determination of the critical temperature with an uncertainty smaller 
than 1 or 2 MeV [16], it is unlikely that these small temperature differences can be directly 
observed. Nevertheless, a trend in the observable signals characteristic of the critical point 
as a function of increasing asymmetry may be detectable. 

The evolution of the system is a complicated process that involves hydrodynamic fiow 
and expansion, together with dissipative effects and nonequilibrium processes like nucleation 
and fragmentation. As the system expands, internal energy is transformed into collective 
motion, which manifests itself as local fiow velocity. If the expansion is highly damped, the 
energy of motion is rapidly transformed back into internal energy; the resulting expansion is 
therefore slow and proceeds along an equipotential surface. The expansion continues until 
p = 0, after which the warm nuclear fiuid remains at rest and evaporates particles until 
it is cool. In contrast, if the motion is undamped, the expansion carries the system past 
hydrostatic equilibrium, where it begins to slow down as the energy of motion is returned 
to internal energy. Since there is no dissipation, the expansion is isentropic, and the motion 
is bounded by the equipotential surface of the initial hot configuration. 

As noted earlier, intranuclear cascade calculations [8,14] suggest that the expansion will 
be isentropic. The system will expand adiabatically into the coexistence region and become 
either superheated (for small specific entropy) or supersaturated (for large specific entropy). 
If the expansion is slow enough to allow for nucleation, bubbles of gas form or the system 
separates into droplets and vapor. If the nucleation is relatively slow, which is more likely, 
and the expansion halts before a spinodal is reached, the direction of motion is reversed and 
the system vibrates, ultimately evaporating neutrons to cool down. In contrast, if the system 
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crosses a spinodal, fragmentation occurs; the relevant mechanism for the decomposition (that 
is, which spinodal is the relevant one) depends on the relative rates for isovector, isoscalar, 
and thermal fluctuations. 

With our simple picture that focuses on equilibrium states, it is impossible to say any- 
thing definitive about the dynamics of the phase separation. This is especially true since it 
is unlikely to occur through a sequence of equilibrium configurations and will involve instead 
superheating, supercooling, or spinodal decomposition. Nevertheless, the dynamics of the 
phase separation is intimately linked to the character of the phase transition; for example, 
the critical exponents are expected to play an important role [6,f6]. Since the phase transi- 
tion in asymmetric matter is not of the usual van der Waals type, there could be observable 
consequences in the signals used to detect the phase separation, such as the distribution of 
emitted mass fragments. In particular, unlike the van der Waals case, the pressure, tem- 
perature, density, and concentration of the gas and liquid phases can change throughout 
the phase separation process. Moreover, since the spinodal decomposition occurs generally 
through a diffusive instability, any microscopic model of the fragmentation process must 
allow the gas and liquid phases to have different concentrations. The extension of existing 
models of nucleation, fragmentation, or percolation (see, for example, Refs. [16,50-52]) to 
allow for this new degree of freedom is an important topic for future study. 

Finally, we discuss how the properties of the phase transition depend on reasonable 
changes in the nuclear matter symmetry energy and compressibility. 
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FIG. 22. Pressure as a function of density at fixed temperature for different values of the 
proton fraction. The nuclear symmetry energy 04 = 30,35,40MeV for the solid, dotted, and 
dashed curves, respectively. 
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FIG. 23. Binodal sections at T = 10 MeV for different 04. Tlie nuclear symmetry energy 
04 = 30, 35, 40 MeV for the solid, dotted, and dashed curves, respectively, and the critical points 
(CP) are indicated in each case. 

The symmetry energy can be varied by changing the coupling , and we chose the values 
04 = 30, 35, and 40 MeV for illustration. (All other couplings in Table I were held fixed.) 
In Figs. 22 and 23, we show the variations in the pressure and in the binodal surface as the 
symmetry energy is changed. As expected, as the symmetry energy increases, the minimum 
y (maximal asymmetry) at which phase separation occurs increases, but only by roughly 
10%. The critical pressure at fixed temperature increases somewhat 40%), but the 

shape of the binodal surface is qualitatively similar; the increased symmetry energy evidently 
has more effect on the gas phase than on the liquid phase. This result is consistent with 
our earlier discussion that it is energetically favorable for nuclear matter to separate into 
a less asymmetric liquid and a more asymmetric gas. As the symmetry energy increases, 
one finds a larger region of configuration space where this phase separation is favorable, 
particularly for retrograde condensation. The critical temperature at j/ = 0.35 decreases 
from r,(0.35) = 13.77 MeV for 04 = 30 MeV, to r,(0.35) = 13.71 MeV for 04 = 35 MeV, to 
rc(0.35) = 13.66 MeV for 04 = 40 MeV, which is a relatively small amount. These modest 
changes as the symmetry energy is varied support the claim made in the Introduction that 
our simple model of the isovector mean-field dynamics is adequate for a discussion of nucleus- 
nucleus collisions. 

We also studied the binodal surface for compressibilities Ky^ = 200, 250, and 300 MeV. 
(We return to our original value of the symmetry energy, 04 = 35 MeV.) The qualitative 
structure of the surface was unchanged, and the most significant feature was a shift in the 
critical temperature of symmetric matter from rc(0.5) = 13.35 MeV at Ky^ = 200 MeV 
to r,(0.5) = 15.75 MeV at Ky^ = 300 MeV. Changes in for other values of 0.3 < 
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y < 0.5 were similar. As noted earlier, these variations of 1 or 2 MeV are unlikely to be 
observable, due to statistical fluctuations near the critical point. We therefore conclude that 
the observation of the liquid-gas phase transition in nuclear collisions is unlikely to provide 
definitive information on the nuclear matter symmetry energy or compressibility, although 
the critical temperature is apparently more sensitive to the latter. 
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VI. SUMMARY 



In this paper we studied the hquid-gas phase transition in warm, low-density nuclear 
matter as a function of the proton fraction. The new ingredient in the analysis was the 
careful treatment of the two conserved charges, baryon number and isospin, which shows 
that the phase transition does not exhibit the usual van der Waals behavior. An examination 
of the stability criteria on the free energy, together with the charge conservation laws and 
Gibbs' criteria for phase equilibrium, reveals that the system can be specified by the same 
number of input variables (temperature and densities) regardless of the number of phases. 
The equilibrium conditions determine the region in parameter space where separation into 
two phases is energetically favorable, as well as the boundary of this region, the binodal 
surface. This surface is two-dimensional for a two-phase system with two conserved charges, 
in contrast to the familiar one-dimensional surface when there is but one conserved charge, 
and this leads to qualitatively new behavior, such as retrograde condensation. For specified 
input variables, the number of equations determining the properties of the phases is equal 
to the number of unknowns, allowing for an unambiguous Maxwell construction of the 
equilibrium state at any point during the transition. In contrast to the usual van der Waals 
case, we learned that in general, the pressure, temperature, density, and concentration of 
both the gas and liquid phase can vary throughout the transition. Moreover, both the 
Gibbs free energy and entropy are continuous throughout the transition, showing that it is 
second-order (by Ehrenfest's definition) rather than first-order. 

To apply these results to nuclear matter, we used a relativistic mean-field model involving 
the interaction of baryons with scalar and vector fields. This model allows for an accurate 
description of the bulk properties of nuclei and of symmetric nuclear matter, which lets 
us calibrate the model and then extrapolate to subnuclear densities and arbitrary proton 
fraction. Although our discussion focused on the qualitatively different aspects of the liquid- 
gas transition in binary systems, the quantitative numerical results are obtained with an 
equation of state that is as accurate as any currently available. By studying reasonable 
variations in the nuclear matter symmetry energy and compressibility, we found that these 
are unlikely to have significant qualitative impact on the signals of the phase transition. 

This thermodynamic mean-field model was then used to study the warm nuclear mat- 
ter that will be produced in energetic heavy-ion collisions. Although the assumption of 
thermodynamic equilibrium oversimplifies the collision dynamics, we believe it is useful for 
providing a concrete description of the bulk properties of the warm matter and for examining 
qualitative features that should have remnants in more microscopic calculations. For exam- 
ple, there are several significant differences between the phase diagram for an asymmetric 
system and that for symmetric matter. First, the critical temperature Tc is different for 
different values of the proton fraction j/, and for j/ 7^ 0.5, Tc is not the maximal temperature 
at which the phase separation can occur. Second, the most extensive spinodal determining 
the instability boundary in an asymmetric system is determined by a diffusive (chemical) 
instability, rather than by a mechanical instability signaled by a negative compressibility or 
by a thermal instability signaled by a negative heat capacity. Thus, in general, the system is 
unstable to isovector modes of separation rather than to isoscalar modes. We also observed 
that in symmetric matter, the isothermal spinodal determines the region of instability for 
either isothermal or adiabatic expansion. (The diffusive and isothermal spinodals coalesce 
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as J/ — )• 0.5.) 

Finally and most importantly, the dimensionality of the phase-separation region is larger 
in an asymmetric binary system, which implies that the phase transition is continuous^ 
and which allows the thermodynamic properties of the phases to change throughout the 
transition. Although the changes in the thermodynamic variables are small for realistically 
observable systems (0.3 < j/ < 0.5) in the model studied here, the resulting spread in the 
thermodynamic variables will increase the variations expected from a study of statistical 
fluctuations alone. Moreover, the increased dimensionality of the phase-separation problem 
could generate significant changes in the observables calculated in a more microscopic treat- 
ment, for example, one based on fiuctuations and spinodal decomposition. At the least, one 
must allow for different proton concentrations in the liquid and vapor when dealing with 
asymmetric systems, rather than including only fiuctuations in the density. Such micro- 
scopic calculations remain to be performed and represent important topics for future work 
on this problem. 

ACKNOWLEDGMENTS 

We are pleased to thank our colleagues G. F. Bertsch, J. F. Dawson, R. J. Furnstahl, 
S. V. Gardner, C. J. Horowitz, and H-B. Tang for useful discussions and comments. This 
work was supported in part by the Deutscher Akademischer Austauschdienst (DAAD) and 
by the U.S. Department of Fnergy under contract No. DF-FG02-87FR40365. 



46 



APPENDIX: THE LOW-DENSITY LIMIT 



Here we derive some useful results for computing the low-density limit of the equation of 
state defined by Eqs. (30) and (47). For densities p ^ 0.005 fm~"^, the results in Sec. Ill lead 
to numerical inaccuracies that can be overcome by making an explicit low-density expansion. 
In general, one has to distinguish between two limiting cases: (z) — )• at finite T and (zz) 
the nonrelativistic limit. We will consider these two cases in turn. 

At fixed temperature, zero density requires Vp^n 0. We therefore consider 

= ez>p,„ and M* = M* + e"AM* for e ^ , 

where Mq is the effective mass at zero density. Using the recursion relations in Eqs. (39)- 
(42), expansions of the densities are readily obtained: 



P = :^{i^P + ^n) [2Hs{0, M*) + M*^i7i(0, M*)j + 0{e') , (AI) 
e 



= —2i^v - ^n) {2HsiO,M*) + M*'HriO,M*)) + Oie') , (A2) 
M* 

2^Hs{0,M;) + O{e') , (A3) 



Ps 



7r2 



which gives 



y = -(l + ^^^^ ) for ^ ^ . (A4) 



2 \ Up -\- 

The effective mass Mq is determined by the corresponding limit of Eq. (43): 

K A M* 

^ (M - M*) + — (M - M; f + 7-AM- M; f = 2^Hs{0, M*) , (A5) 

which is independent of the proton fraction y. By inspecting the lowest-order corrections in 
e, it follows that a = 2, i.e., 

M* = m; + o{p^) . 

Moreover, by observing that 

W=p{p) and R=p{p), 
we conclude that all isotherms with < y < 1 approach a common limit at zero density: 

P = :t\^5(0,M*)-^(M-M*)^--^(M-M*)^ 

^-{M-M;r + 0{p'). (A6) 



2ig 

At moderate temperatures [T ^ Af/IO), this vacuum pressure is very small. 

We turn now to the nonrelativistic limit, which is encountered in particular in adiabatic 
processes. This limit can be characterized by 

— )• and M* /T — > oo , but rjp^n = {i^p,n — M*)/T finite. 
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It is a straightforward task to derive asymptotic expansions for the integrals (37) and (38): 
Gn{p.M) ^=^^ i^'^^'(^F„/,_,(^) + („ + 2)^F„/,(^) + 0(rVM^)) , (A7) 



M/T^oo 

with 



T{2MTtl'-'^ {f^„.M + (n - 2)^F„/2(^) + 0{T' jM^)]^ , (A8) 



/•CO 

Fa{v) = / dx- 



lo 1 + e^-*? 

To lowest order, this yields the familiar Fermi gas results for the densities: 

p^p^^i^^Mp.^ [Fi/2(?7,) + Fi/2(?7.)] (A9) 

and for the entropy 



(7 



^ {5[F3/2{Vp) + F3/2{Vn)] - 3[v,F,/2{Vp) + VnF^f2{Vn)]) • (MO) 



The effective mass in these expressions has been replaced by the nucleon mass Af , since 



M*^M-^p, (All) 
mj 

which follows from the relation (A9). In an adiabatic process, the constants rjp^n can be 
obtained for a given entropy /baryon s = a / p and proton fraction y using 

a 5 \F3/2irip) + i^3/2(??n)l - 3 [?7pFi/2(?7p) + ?7„Fi/2(??n, , 

p (^^^) 

^ Fi/2{Vp) + Fl/2{Vr, 

2y-I = ^^4^^^^. (AI3) 
By including the lowest-order contributions in the density, the pressure then takes the form 



3M ' {F^/2{np) + F^i2{r}n)fl^ 

vf-.p" - f-.p" + ^(2j/ - Ifp' . (AI4) 
Zmfj Zmj Snip 



As expected, at very low densities the pressure is dominated by the ideal-gas term with 
adiabatic index 5/3. The nonlinear couplings give rise to higher-order contributions in the 
density and can be neglected. 
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